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Abstract 

Small angle neutron diffraction experiments are analyzed using recently de- 
veloped and properly generalized one-field effective free energy method. In 
the case of experiment of Keimer et al on YBCO, we show that the fourfold 
symmetry of the underlying crystal is explicitly broken, but the reflection 
with respect to the [110] and [110] axes remains a symmetry. The vortex 
lattice also becomes generally oblique instead of rectangular body centered. 
Unexpectedly rich phase diagram is described. 
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There are growing evidences that superconductivity in layered high Tc cuprates is largely 
due to the (i(a;2_y2) pairing with small mixing of s-wave component The uncon- 

ventional pairing mechanism makes an impact on the single vortex and the vortex lattice 
structure. Recent studies on the detailed structure of the Abrikosov vortex lattice in YBCO, 
using small angle neutron diffraction and tunneling spectroscopy |]^, show clear devi- 
ations from the standard triangular lattice. It is natural to try to explain these deviations 
theoretically with modified phenomenological Ginzburg - Landau (GL) theory. To investi- 
gate d-wave superconductors with s-wave mixing, Ren et al. P] and Soininen et al. both 
derived an effective GL type theory using two order parameters: s and d. From this effective 



action, or more fundamental equations |]T0|, one obtains a characteristic four-lobe structure 
for an isolated vortex and its associated magnetic field The fourfold vortex core struc- 
ture comes into conflict with the high symmetry of the triangular lattice and can distort it 
at already accessible fields much lower than Hc2- The vortex lattices obtained within this 
approach are basically centered rectangular lattice with chains of vortices oriented along 
crystalline axes [100] and [010] (see Fig. 1). They spontaneously break the fourfold ro- 
tational symmetry (i.e., two different lattices related by 90° rotation), but preserve the 
reflections with respect to the axes [100] and [010]. 

These predictions come close to results of some experiments [§|,05 but clearly disagree 
with those of [Q. According to the interpretation given in the centered rectangular vortex 
lattice gets rotated by 45° with respect to the crystalline axes (see Fig. 2a), i.e., the chains 
of vortices lie along the diagonal directions [110] and [110] instead. A recent theoretical 
study by Ren, et al. [|l^ has considered explicit breaking of the fourfold symmetry within 
the two field framework. Their results however remains qualitatively the same as the case 
with fourfold symmetry - only centered rectangular nonrotated vortex lattices are obtained. 
So far, there is no theoretical interpretation for the lattice data observed in ||^. We shall 
provide such an interpretation in this letter. Our answer is different from that provided 
in 1^, however, the results can still be derived from the GL theory with proper fourfold 
symmetry breaking terms. 

In this work we adopt a recently developed one field effective theory, first introduced 
by Affleck et al |]13[ in which they work mainly in the London limit, and later by us |T^ 
for static and moving vortex lattices near ifc2- Most of the above mentioned results can be 
reproduced in this much simpler formulation in which only the field d is introduced and the 
theory is based on the following D^h symmetric free energy 

Feff[d] = - + m' - vd* (n^ - Ulfd, (1) 

where 11 = —iV — e*A. The last term which we call F^^ parametrizes the breaking of 
full rotational symmetry down to D^h and can be treated as a perturbation. Near Hc2, 
the linearized equation in the one field approach can be solved perturbatively in rj, which 
allows one to easily generalize the description of the centered rectangular lattices to the 



most general oblique lattices [|14[. This will be crucial in the present work in which these 
more general lattices are indeed the ground state in some cases. 

Note that the contributions to the coefficient 1] might not only come from the d-s mixing 
which always gives a positive r], but also from other sources []I1[. The possibility of having 
negative rj will be discussed later. It is also important to realize that since this formula- 
tion utilizes only the symmetry properties, it can be applied to the conventional type II 



2 



superconductors with symmetry as well. In this case, t] is proportional to the angular 
average of products of Fermi velocities on the Fermi surface, describing the deviation of the 
Fermi surface from a perfect sphere |^ . The effective strength of F4^d can be characterized 
by a dimensionless parameter rj' = rjiride*!! Using the free energy in Eq.(|I]), one finds 
centered rectangular vortex lattices (see Fig. 1) with the angle a directly related to the 
coefficient rj' . The lattice becomes square when rj' exceeds a critical value 77^ = .0235 [jl4 . 
This can accommodate the tunneling spectroscopy data of [0 and the SANS data of ||^ for 
YBCO, as well as a recent decoration and neutron scattering data for a low T^. material 
ErNi2B2C |]T3[. The analysis presented in |]T^ indicates that the precise SANS data of ||T^ 
unambiguously shows that for large the vortex lattice becomes a square one, exhibiting 
perfect D^^h symmetry. The less precise data of gives an angle a ~ 54°, which corre- 
sponds to rj' 0.019 and is in the centered rectangular phase. These two experiments on 
two different samples both seem to show manifestation of D^^ symmetric GL free energy and 
correspond to its two different phases. The transition from the centered rectangular vortex 
lattice to the square lattice was observed in ErNi2B2C and has not been observed yet 
in high materials. 

In order to explain the data in , we now generalize the formalism to include terms which 
break the D^^ symmetry. This can be also motivated by noting that in many high cuprates 
the Di^h symmetry is not exact. For example, the CuO chains in YBCO breaks the fourfold 



symmetry down to twofold |]T2| . Up to (scaling) dimension three, there are two possible terms 
that break fourfold symmetry: Fx2_y2 = —fid*{Ily — Ill)d and F^y = —Xd*{]lxlly + ]lyllx)d. 
The first term Fx2_y2 describes the asymmetry between [100] and [010] axes and has the 
refiection symmetries x — ^ 
already been considered in 
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This term has 



121. The second term F, 
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on the other hand, preserves the 
refiection symmetry with respect to the [110] and [110] directions, that is, x — ^ y, ?/ — > x 
and X —>■ —y,y —>■ —x. In the BCS theory, the presence of the second term requires that 
the shape of the Fermi surface also breaks the and ay symmetries . Since this is quite 
unlikely, we do not expect that it will occur in the conventional superconductors. We will 
find, however, in the case of Keimer et al's SANS experiment, the F^y term is required to 
explain the data. 

The method of calculation is quite analogous to that of the 77 correction explained in 
I^ , so here we just present the result. Let a,b be the two lattice constants and a be the 
angle between the two basis vectors (Fig.l). It will be convenient to introduce the complex 
variable ( = -e*" = p + ia. The angle between the vortex lattice and the crystalline lattice 



will be denoted by ip. The Abrikosov's Pa = (MI"^) / (MI^) is then given by: 



(3a (p, (t) = (p, a) + ^ Re 
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where (p, a) can be found in standard textbook or in [^ . All the three anisotropic 
corrections are collected in the prefactor: 

G{n) = r]'e^''^{6An^a^n^ - ASTran"^ + 3) 
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+4/iV^^(87r(Tn2 - 1) 

+4A'e2*(^+'^/^)(8W-l) (3) 

where /x' = /im^, A' = Am^. 

The term Fy2_^2 in the effective energy preserves the symmetries of the centered rectan- 
gular lattice and is therefore not expected to produce interesting qualitative effects, so we 
will drop the fi' term in the following discussions, however it is understood that in making 
quantitative comparison with data, the n' term may have to be included. The remaining 
correction to Pa summarized in G{n) has two parts: the first one comes from the fourfold 
symmetric term F^d and has e^**^ angular dependence. The second one has e'^^^'^^'^^^^ angular 
dependence and comes from F^y. It is this conflict between the two contributions that gives 
rise to the observed diffraction pattern. Either F^d or Fxy alone will give reflection invariant 
lattices, i.e., rectangular body centered lattice aligned along [100] or [110], respectively. The 
lattice structure is determined by minimizing (3^ with respect to p, a and if numerically. 
One obtains generally nonrectangular oblique vortex lattices. It differs markedly from the 
D^h symmetric case. 

Fig. 2(b) shows the diffraction pattern and the corresponding lattice structure that we 
obtained at rj' = 0.019, A' = 0.04. (Note that in 2D the reciprocal lattice is nothing but a 
rotation of 90° of the real lattice). In one of the diffraction patterns. Fig. 2 (a), one sees clearly 
two large peaks in the [110] direction and four weaker points on both sides of the [110] line, 
giving totally 10 points. This was interpreted in Ref. §] as a nearly rectangular lattice with 
one of the basis vector lying on [110], together with its reflected version. (Presumably, the 
two lattice orientations are degenerate ground state and show up simultaneously as different 
domain in the sample). The points on the [110] line then coincide and produce constructive 
interference. In comparison, in the previous calculations, because the reflection symmetries 
ax and ay are preserved, one always obtained rectangular lattices which are aligned to either 
[100] or [010]. They possess twofold symmetry and upon reflection one is unable to produce 
different lattices. As a result, there will be only 6 points on the diffraction pattern and can 
not account for the data. 

We notice an important difference here: The off diagonal points (the four weaker points) 
are not really on the line parallel to [110] as Keimer et al. claimed. One might hope to 
tune the parameters such that when the two points on [110] merge into one, the off diagonal 
points will align themselves as well, but this is not the case. In fact, they will also merge with 
each other, and there will be no splitting anymore. If one look carefully at their contour plot 
it is possible to tell the difference. Furthermore, the lattice we obtained is not rectangular, 
this is consistent with their possibly 5% difference between the length of the two primitive 
basis vectors. 

The vortex lattice phase diagram in the [i]', A') plane is presented in Fig. 3. Since 
changing the sign of A' only reverses the roles of [110] and [110] axes, it suffices to show only 
the positive A'. First, consider the D^^ symmetric case with A' = 0. Then rj' = corresponds 
to the conventional triangular lattice with no special orientations. For t]' < rj'^ the lattice is 
centered rectangular aligned to [100] and [010] with double degeneracy (related by reflection 
about [110]). Increasing rj' elongates lattices along either [100] or [010] so that when rj' > r]'^, 
the two degenerate lattices both becomes square and the full symmetry is restored. 

For A' > 0, there are three phases and two phase transition lines. The lattice, compared to 
the corresponding A' = case, can in general be thought of as resulting from a deformation 
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in the [110] direction. (For A' < 0, the corresponding deformation will be in the [110] 
direction). The lattice is centered rectangular for smaller values of , while it is rectangular 
(not centered) for larger values of 77'. Symmetry of the unique ground state in each of these 
two cases is larger than that of the free energy. There is, however, no direct phase transition 
between them. Instead, in the region between these two phases bounded by A'i(r7') and ^^k^^i 
there is a less symmetric phase in which ground states are doubly degenerate. This comes 
from a nontrivial competition between and F^y The two degenerate lattices are also 
related by the reflection about [110] and are generally oblique. We see that the data in [^ can 
be fitted into this phase. The transition line A'^(r7') starts from the origin and monotonically 
increases with r/', while A2(?7') starts from (^7^,0) and also increases monotonically. A2(?7') 
appears to approach A'i(?7') asymptotically. Since r]' is proportional to the magnetic field if, 
one immediate implication of this phase diagram is that, for a given sample, by increasing 
H one should encounter two phase transitions. This prediction can be tested directly by a 
number of experimental techniques. 

We would like to briefly describe here another, rather exotic possibility. The one field 
approach allows one to consider the negative r/ case. This cannot be obtained from the two 
field formulation in which 77 is always positive if we only assumes one critical temperature 
I^J^. However the possibility of negative // cannot be ruled out theoretically. In the one 
component theory with exact fourfold symmetry, the negative r] is equivalent to the d^y 
pairing, while in the BCS theory, it could happen if the Fermi surface is elongated along 
y = ±x direction. When t]' is negative, the minus sign replaces (p in the exp {4:i(f) factor in 
Eq. (H) by (/? ±45°, and then both and F^y will prefer the diagonal direction. As a result 
there will be no competition and we will always get rectangular body centered lattices along 
[110] or [110]. 

In conclusion, we have investigated the effects of explicit fourfold symmetry breaking on 
the Abrikosov lattice structure with the one field formulation . The complete phase diagram 
was constructed. We found quite rich phase diagram with three different phases separated 
by two phase transition lines. The vortex lattice observed in Keimer et al.'s experiment ^ 
can be accommodated in the new phase diagram . It turns out that the vortex lattices are 
no longer centered rectangular, but rather general oblique ones. Other experiments fit quite 
well into the D4/J symmetric phase in which the triangular to square phase transition takes 
place. 

Authors are very grateful to CC Chi, M.K. Wu, W. Yang, and S.Y. Hsu. This work 
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Figure Caption: 

Fig. 1 The body centered rectangular lattice obtained in the fourfold symmetric case, the 
two lattices (a) and (b) are related by a rotation of 90" or reflection about the [110] axis. 

Fig. 2 Keimer et al's SANS diffraction pattern and two different interpretations, (a) Keimer 
et al's interpretation, and (b) The interpretation given in this paper. 

Fig. 3 The phase diagram for the vortex lattice structure as a function of the four fold 
anisotropy parameter vj and the two fold anisotropy parameter A' . 
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